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Abstract 



We study the exclusive decay of J/^' into a lepton pair combined with light 
Ph. hadrons in the kinematic region, specified by that the total energy of the light 

D . hadrons is much smaller than rUc, the mass of the c-quark. In this region, the 



nonperturbaive effect related to J/^ and that related to the light hadrons 
can be separated, the former is represented by a NRQCD matrix element, 
while the later is represented by a matrix element of a correlator of electric 
chromofields. The results are obtained in a axial gauge by assumming that 
contributions from two-gluon emmssion are dominant. But we can show that 
these results can be obtained without the assumption in arbitrary gauges. A 
discussion of the results are presented. 
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1. Introduction 

Several ten millions J/\& events will be collected with Bejing Spectrometer (BES), this 
provides many opportunities to study different decay modes of J/'^ with high statistics. 
In this work we propose to study the exclusive decay into a lepton pair combined light 
hadrons with a total energy much smaller than the mass of J/^ ^ this decay will be studied 
experimentally at BES [|^ 

Many interesting properties of QCD can be learned from studies of J/\l/ or a quarkonium 
system. J mainly consists of a c-quark and its antiquark and those quarks move with a 
small velocity Vc inside J/\E' in its rest-frame. Hence an expansion in Vc can be employed 
to describe properties of J/\E' . The expansion can be systematically performed in the 
framework of nonrelativistic QCD(NRQCD) [0] for inclusive decays and for production rates. 
In this framework the inclusive decay of J/\l/ can be imagined at the leading order of Vc as 
the following: The c- and c- quark in J/\l/ has certain probability to be freed at the same 
space point and this cc pair decays subsequently. The probability has a nonperturbative 
nature and is at order of v °, while the decay of the cc pair can be treated with perturbative 
QCD. Various effects at higher order of can be taken into account, e.g., relativistic effect, 
and the effect of that the freed cc pair does not possess the same quantum numbers as those 
of J/^>. 

For the exclusive decay we will study, this interpretation still can be given at the leading 
order of Vc- Then the decay is mainly of the decay of the cc pair and it can be thought that 
the pair emit some soft gluons and annihilates into a virtual photon. The soft gluons will be 
transmitted into the light hadrons with a total energy which is small. This interpretation 
implies that the nonperturbative effect related to J /'^ and that related to the light hadrons 
can be separated, a factorization of the S'-matrix element for the decay can be performed. 
It should be noted that the factorization here is performed at the amplitude level for the 
exclusive decay. It is not like the factorization for inclusive decays proposed in p. We will 
assume that the factorization can be performed and then try to obtain a factorized form for 
the decay amplitude. A proof for the factorization in our case may be done as done for other 
exclusive process Q. In our case one may also take the J/\l/ as a nonrelativistic bound-state 
of a cc pair and use a wave-function for the bound-state. Then the wave-function at the 
origin can be written as NRQCD matrix elements defined in and our results can be 
obtained. 

For the light hadrons we will first assume that the soft gluons consists only of two gluons 
and take the gauge = 0, where G^{x) is the gluon field. Because the soft gluons 

have smaller momenta compared with the c-quark mass mc, an expansion in can be 
used. At the leading order, the S-matrix element for the decay is suppressed by m~^, this 
may be understood in the framework of HQET 0. In the axial gauge the leading term of 
HQET does not contain gluons, hence an emission of one or two gluons is suppressed by 
m~^. At the leading order of the decay amplitude can be factorized into three parts: 
The first part is a NRQCD matrix element representing the bound-state effect of J/\l/, the 
second part is a matrix element of a correlator of electric chromofields, which characterize 
the soft-gluon transition into the light hadrons. The third part consists of some coefficient. 
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However, the assumption of two-gluon-emission is not justified, because it sounds that we 
use perturbative theory for soft gluons. But we can show that without the assumption one 
can also obtain the same results only by using the expansion in m~^. Hence the results are 
nonperturbative and the number of emitted gluons is not restricted in an arbitrary gauge. 
In our work we only consider the emission of soft gluons which are transmitted into the 
light hadrons. Effects of gluons exchanged between quarks, between gluons and between 
quarks and gluons are not considered. If the exchanged gluons are hard, then their effect 
can be calculated with perturbative theory and it results in a correction to our results at 
higher orders of a^. If the exchanged gluons are soft, their effect is nonperturbative and 
may be factorized into nonperturbative matrix elements provided that the factorization can 
be proved. It is beyond this work to consider the effect of exchanging soft gluons. 

Typically, the light hadrons are two pions. Hence the decay studied here has similarities 
comparing with the decay J/\E' + vr + tt, which is studied in where the soft-gluon 

transition can be described by matrix elements of local fields in the decay width. In our 
case the transition is characterized by a matrix element of nonlocal fields. However, they 
are related and we will show this in detail. 

Our work is organized as the following: In Sect. 2 we derive our results for the decay 
with the assumption of the two-gluon emission in the axial gauge and discuss the relation 
to the decay \E'' ^ J/\E' + vr + tt. In Sect. 3 we will derive the results in arbitrary gauge and 
without the assumption, i.e., the number of emitted gluons is not restricted. In Sect. 4 we 
summarize our work and discuss our approach. 

Throughout of our work we take nonrelativistic normalization for the J/^P state and for 
c- quark. 

2. The Decay with Two-GIuon Emission 

We consider the exclusive decay of J/\l/ in its rest-frame: 

J/^(P) ^ i+{pi) + r{p2) + light hadrons, (1) 

where the momenta are given in the brackets and the light hadrons are specified hadrons 
with the total momentum k, whose component is much smaller than m^. Typically, they 
are two soft pions. At the leading order of QED the S'-matrix element for the decay is 

{f\S\i) = -le^QcL, ■ 4 / dhe''i-'{LH\-c{z)rc{z)\J/^), (2) 

q J 

where Qc is the electric charge of c-quark in unit of e, c{x) is the Dirac field for c-quark, 
{LH\ stands for the final state of the light hadrons and 

q = Pi+P2, 

Lf, = u{p2h^,v{pl). (3) 

In Eq.(2) we only take the part of c-quark in the electric current into account. Assuming 
that only two gluons are emitted by the c- or c-quark, we obtain 



3 



{LH\T [c{x)j ■ G{x)c{x)c{y)j ■ G{y)c{y)c{z)rc{z)] \J/^), (4) 

where G{x) is the gluon field. Using Wick-theorem we can calculate the T-ordered product 
and we only keep those terms in which one c-field and one c-field remain uncontracted. Then 
the matrix element takes a complicated form and can be written in a short notation: 



{f\S\t) = -i-e^Qc9lL, -^j d'xd'yd^zd'xid'yie'''-^{LH\Gl{x)Gl{y)\0) 

(0|c,(xi)Q(t/i)|J/*) •M;7''^'(a;,y,xi,yi,^), (5) 



where MjJ'^'^'"^(a;, y, xi,yi, z) is a known function, i and j stand for Dirac- and color indices, 
a and h is the color of gluon field. The above equation can be generalized to emission of 
arbitrary number of soft gluons, then the S'-matrix element is the sum of the contribu- 
tions with 2,3, ■ ■ ■ soft gluons and in each contribution there is the same matrix element 
{Q\cj{x)ci{y)\J 1'^). For this matrix element the expansion in Vc can be now performed, the 
result is: 

(0|c,(x)Q(y)| J/*) = -i(P^yp_),,(0|xt<7V| J/*)e-^^-(^+^) + 0{vll (6) 

where x^('0) is the NRQCD field for c(c) quark and 

P± = (1 ± r)/2, 

p'^^K, 0,0,0). (7) 

The leading order of the matrix element is 0{v^), we will neglect the contribution from 
higher orders and the momentum of J/^' is then approximated by 2p. It should be noted 
that effects at higher order of Vc can be added with the expansion in Eq.(6). Taking the 
result in Eq.(6) we can write the ^'-matrix element as: 

{m^) = -^2^e^ge5.'(27^) V(2p - k - q)L, ■ 1 

(0|xVV|J/*)2Y d'x j ^^e^^'^^-^{LH\Gl{x)Gt{-x)\0)R^'''^'{p,k,q,). (8) 

To obtain the equation we have used the color-symmetry and the translational covariance. 
The quantity R takes the form 

1 1 

7 • {-P + 1^ ~ ?i) -mc + iO+ ^ ■ (k-p) -mc + iO+ 

7 • {-P + - Qi) - + ?0+ 7 • (p - |/c - gi) - mc + iO+ 

u,u ^ li,qi^ -qi)} (9) 
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The physical interpretation of R is that it is the amphtude for a ^Si cc pair emitting two 
soft gluons and a virtual photon, and c and c have the same momentum p. We take the 
gauge 

G^x) = 0, (10) 

and in this gauge the electric chromofield is given by: 

E{x) = E"(x)T" = doG{x). (11) 

The quantity R depends on p, k and qi. In the kinematic region we consider, A; is a small 
vector. The dominant region for integration is characterized by k and by Aqcd, because 
the x-dependence of the matrix clement {LH\ • ■ • |0) is characterized by k and by Aqqd. 
Therefore one can expand R{p, k, qi) in m~^. Keeping only the leading order we obtain 



0\2\ 



{LH\G\x) ■ G"(-x)|0) + O(^) (12) 



ml' 



Substituting the above expression into the ^'-matrix element, the integration over qi and 
over X can be performed. Then we use the following equations 



d''xe^"'-%k + 2q)p{LH\Gl{x)Gl{-x)\Q) = -2i j d^xe^"' '' {LH\Gl{x)dpGl{-x)\0) 
j (iW^^-^(A; - 2q)p{LH\Gl{x)Gl{-x)\^) = -2i j d'^xe^''i-^LH\dpGl{x)Gt{-x)\0), (13) 



and we obtain 



n'^(0^\'^^'^(Or. -h-n\T. . ^/n\.A^^o ^ ^ 



{f\S\i) = i-e'Q^g',{2nm2p - k - q)Lp • ^(0|xVV| J/*) • — ^^o)2 

. / ^ 1 \ j dte'^^'HLH\W{t,Q) -Wi-tMQ) 

i 27rl + r-i0+ 1-r-iO+y ^ \ \^ ) \ ^ )\ / 

+0{^) + 0{v% (14) 

where r is related to q\ by 2q\ — rk^. The term with gf^ is expected in the heavy quark 
limit. In this limit gluons will not change the spin of c- or c quark. The matrix element 
{0\x^a^ip\J/^) is proportional to the spin of J/^, hence the spin of J/^ is transferred to 
the virtual photon. The result is obtained in the axial gauge defined in Eq.(lO), to make 
the S'-matrix element gauge invariant we must add a gauge link between the two operators 
of the electric chromofield. We define a distribution amplitude for the gluon conversion into 
the hght hadrons: 

Ht) = dte'^'°'{LH\B'^{t,0) ■ [PeM-i9s£/x'G'^^{x',0)T'^}l,E\-t,0)\0), (15) 
where P means path-ordering and r*^ is the generator of SU (3) in adjoint representation 
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{T'')ab = -ifabc- 



(16) 



The function h{T) is gauge invariant, it represents the transition of soft gluons into the hght 
hadrons. With the gauge hnk the transition is of two space-hke gluon plus many time- 
hke gluons. It also depends on each individual momentum of each hadron, we denote the 
dependence as /i(r, k). Because of the energy conservation /i(r, /c) = if |r| > 1. With this 
the S-matrix element can be written: 

9 nP^ 1 ^ 



{f\S\z) = ^-eU^7f (2vr) V(2p - k - q)L, ■ ^ {0\x^a'^\J/^) ■ — 



TLHik) = I dr i h(T,k). (17) 

Eq.(15) and Eq.(17) are our main results, derived in the axial gauge with the assumption of 
two-gluon emission, then we add a gauge link to make the S*- matrix element gauge invariant. 
This may be unsatisfied, because the assumption of two-gluon emission sounds that we 
performed an expansion in Qs for soft-gluons and we add the gauge link by hand. It is 
possible that the cc pair emits soft gluons which are all time-like and form the light hadrons, 
and this emission is not suppressed by m~^. This type of contributions is excluded with the 
gauge. Therefore, it is important to have the results derived in an arbitrary gauge without 
the assumption. In the next section we will show that the results can be obtained without 
the assumption and the gauge-link is automatically supplied. 

Before ending this section we would like to discuss the relation between the decay studied 
here and the decay \E'' ^ J/\l' + tt + tt. If we expand the integrand in g° in Eq.(12) and then 
perform the integration we obtain: 

{f\S\i) = 4eU^?,^(27r) V(2p - k - q)L, ■ ^(OlxVVl^/^) 



E 7^JLH\0,^\^) + O(-l) + 0{v% (18) 

n=0 I 

where are well known twist-2 operators 

02n = z'"G'"'°'^(ao)'"G"^'°''^(7M. (19) 

in the axial gauge. With this we can also obtain a series for TiHik): 

TLHik) = 'l.Y.^^{LH\O2^\0). (20) 

It should be noted that T^h is defined at a renormalization scale fi, the evolution of the 
operators with /i is well known, hence one can also obtain the evolution of T^h- The 
same operators appear in the decay \E'' — > J/\E' + vr + vr 1^-0, the S'-matrix element for the 
decay can be written: 
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(J/* + LH\S\^') = Yl C2n{LH\O2n\0), (21) 
n=0 

where we denote the state of the two pions as {LH\. In this decay the coefficient c„ is 
proportional to m^". Hence one can truncate the series as an approximation. In our case, 
the corresponding coefficients are not suppressed by the power of m~^, we can only sum 
the series into a function Tlh- This also prevents us from giving numerical predictions by 
knowing several leading terms in the series. 

3. An exact derivation of the results 

In this section we work in an arbitrary gauge and do not assume of two-gluon emission. 
The number of emitted soft-gluons is not restricted. With the discussion after Eq.(5) in the 
last section and after performing the expansion in Vc, the problem can be formulated as a 
transition of the cc pair into a virtual photon and the light hadrons. The ^'-matrix element 
can be written as 

{f\S\i) = -ie'Q^L, ■ ll(0|xVV| J/*) 

J2 v{p, S2h'u{p, s,) [ d'ze'^-^LH\ciz)rcizMp, s,), c(p, s^)) + 0{vl) 

= -ie'Q,{2nY6\2p - k - q)L, ■ li(0|xVV| J/^) 

o 

^ v{p, S2Wu{p, s,){LH\c{0)Ycmc{P, si), c(p, S2)) + 0{vl), (22) 

S1S2 

where the c- and c-quark has same color and the summation over the color is implied, Si 
or S2 is the spin of the c-quark or of the c-quark, respectively. Now we need to calculate 
the matrix element {LH\c{0)^^c{0)\c{p, Si), c(p, 52))- For this we use LSZ reduction formula 
and obtain: 

{LH\c{0)Ycmc{p, si), c{p, S2)) = -\j d*xd^ye-'P<^+Mp, S2)(^7 ■ - m,) 

■{LH\T {c{x)c{0)rc{0)c{y)} \0){i^ ■ "dy + mMP, (23) 

where Z is the renormalization constant for the field c{x). We can evaluate the matrix 
element with help of the QCD path-integral, and perform first the integration over c-quark 
field we obtain: 

{LH\c{0)rc{0)\c{p, si), c{p, S2)) = ^J d^xd^ye-'P-^^+y^vip, S2)(^7 ■ - m,) 

{LH\S{x, 0)7^^5(0, y)\0) {zj ■ ty + mMP, (24) 

In the above equation the average over gluon fields, i.e., the integration over gluon fields, 
will be taken later, and S{x, y) is the c-quark propagator in a background of gluon fields: 

-^S{x,y) = (S)\T{c{x)-c{y)}\Q). (25) 
7 



We can define wave-function ipcix) and ipc{x) for the c- and c-quark respectively in a back- 
ground of gluon fields: 

ijj^ix) = J d'^ye-'f-y [S{x,y){i-f ■ ^ y + m^)] u{p, Si), 

Mx) = J dSe^^' MP, S2){^l ■ dy - m,)S{y, x) (26) 

These wave functions satisfy the Dirac equations 

{i'y ■ D - mc)il!c{x) = 0, 

V'c(x)(^7-S + ^c) = 0, (27) 

and the boundary conditions 

^ipcix) e~^^'^u{p, Si) for x° — oo, 

ifjsix) e-'P-^v{p, S2) for x^ -00, (28) 

where = + igsG^{x), G^{x) is the background field of gluons and ipD = {D^ip)'^'^^. 
It should be noted that the propagator in Eq.(25) is a Feymann propagator, wave-functions 
defined in Eq.(26) do satisfy the Dirac equation with a background field of gluons, but 
they do not have a simple boundary condition like in Eq.(28). However, it was shown |Q] 
that the Feymann propagator in Eq.(26) can be replaced with a retarded propagator, if the 
background field varies enough slowly with the space-time. Then one can have solutions 
of Eq.(27) and they satisfy the boundary conditions in Eq.(28). In our case we will make 
the expansion in and in this expansion the c-quark and c quark are decoupled, and it 
results in that the Feymann propagator will automatically be a retarted propagator for c 
and for c. Therefore one can obtain a solution from the Dirac equation and the solution 
can satisfy the boundary condition in Eq.(28). It should be also noted that the boundary 
condition in Eq.(28) implies G^{x) —>■ for x° —00. 
With the wave-functions we have 

{LH\c{0)rc{0Mp,si),c{p,S2)) = ^{LHlUOh'Mm) (29) 

and it is gauge- invariant. If we take c-quark as a heavy quark, we can expand the wave- 
functions in m~^. The expansion is the same as that in HQET 0, we write the wave- 
functions as: 

M^) = e-'P-^ \h,{x) + ^7 ■ DtK + 0(m;2)| , 

^.-(x) = e-'P-^ \h,{x) - l^M^h ■ "Dt + 0{m-^)^ , (30) 
where Dj- = (0,D). For the functions hc{x) and hc{x) we have 

{^Do - ^(7 ■ DTf)hc{x) = + 0{m-'), 

Mx){tBo - ^(7 ■ ^T?) = + 0{m-'), (31) 
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and the boundary condition for these functions can be read from Eq.(28). Giving a back- 
ground field of gluons we solve the above equations by expanding these functions in m~^. At 
the leading order c-quark will only interact with time- like gluons and the solution is simple. 
We define 

dtG^{t,x)} (32) 

-oo 

and with the solution at the leading order the wave-functions are 

^^(x) = e-'P-''V{x)u{p) + 0{m-^), 

^g(a;) = e-'P-''v{p)V^{x) + 0{m-^). (33) 
Because V^{x)V{x) = 1 we obtain after averaging the gluon fields: 

{LH\ciO)rcmc{p,s^),c{p,S2)) =0 + O{m;'). (34) 

That means that the leading order of the S-matrix element is at and only time-like 
gluons will not lead to a contribution to the S'-matrix element at order of m°. This is 
in consistency with the results in the last section, where the axial gauge and perturbative 
theory is used. Here we derive this and the results at order of mj^ in an arbitrary gauge 
and without using perturbative theory. 

To calculate the matrix element at order of m~^, we first make a gauge transformation: 

ij'^ix) = V\x)Mx), ^',{x) = Mx)V{x), 
K{x) = V^x)K{x), h',{x) = h^,)V{x), 

G'^ix) = V\x)G,{x)V{x) - -V^{x)d,V{x). (35) 

9s 

After the transformation we have 

{LH\c{0)rcmc(p,si),c(P,S2)) = ^{LHmO)ri^c(m. (36) 

With this transformation one can show G'^{x) = 0. The wave functions and ip'^ have 
the same expansion in as those in Eq.(30), with the replacement that all functions 
are replaced by functions with a prime. The functions h'^{x) and h'^{x) satisfy the same 
equations in Eq.(31), where the covariant derivative is 0^ = 8^ + igsG'^{x). Solving these 
equations with the boundary condition we have 

K{x) = {l - ^ f ^ ^^(7 • DT{t, x))^| u{p) + 0{m:% 

h',{x) = v{p) |l - ^ £^ dt{^ . Drit, x))^| + Oim:'). (37) 

With these solutions we have 
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4^(a;) = e-^^-^v{p) |l - ^ £^ ^^*G'(t, x) • G'(t, x) + • • - j + OK'^). (38) 



In the above equations • • • in the bracket {} will not lead to contributions at the order we 
consider. Using these results we have 



B'^ = E v{p,S2Wu{p,s^){LH\^-Ml'Mm) 

Sl,S2 

= ^^^(Li/|G'"(0)G'"(0)|0) r dte''''' + 0{mf), (39) 

where we have averaged the gluon fields. It should be noted that is the total energy of 
the light hadrons, which forms an asymptotic final state, hence should be understood 
as A;" — iO"^. With this in mind the integration over t can be performed. One can also use 
^-function 

^ i ^ = (40) 

to perform the integral. For the transformed gauge fields G'^ we have 

E'(x) = E'"(x)r« = aoG'(a;), (41) 
then the gluon field G' can be related to E' by 

d^E'(e,x)=/ ^(i-OE'(e,x). (42) 

-00 J —00 

Using this relation and translational covariance we have 

/OO 
dtdTe'^^e{-t -T)9{t- T){LH IE'" {t)E"'{-t)\0) 
-00 

= 2 r dt%^-^ ' ' 



00 27r oo-lk^ + m (jj + ^k^-iO+ 



■{LH\E"'{t)E'^{-t)\0). (43) 
Changing the variable 2u! — t and we obtain: 



= i-e^ge^,^(27r)^5^(2p - A; - q)L, ■ ^l-(0|xVV| J/*) 



+0(-^) + 0(vl), (44) 
mi 
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Expressing E' with E we realize that the gauge hnk In Eq.(15) is automatically generated, 
{LH\E'%t,0) ■ E"^(-t,O)|0) = 

(L/7|E"(t,0) ■ [Pexp{-z(7,|_*^rfx°G°'^(x°,0)r^}],,E''(-t,0)|0) (45) 

while in Eq.(15) and Eq.(17) we sandwich it by hand. 

The last factor which needs to be considered is the renormalization constant Z for wave- 
function. Following PI we consider the current 

J^ix) = c{x)Yc{x). (46) 

The current J^{x) is a conserved current and it will be not renormalized. Hence we have: 

{c{p)\J^mc{p))=u{p)l'u{p). (47) 

Using the same techniques in this section, one can derive 

Z =l + 0{m-^). (48) 

With the results here we find that the S'-matrix element obtained here is exactly the same 
obtained in the last section. We emphasize here that it is derived here in arbitrary gauge 
and without the assumption of two-gluon emission. Hence the results are nonperturbative. 

With the given S'-matrix element, one can calculate the differential decay width. We 
obtain for unpolarized J/\l/: 

^ =«Q^^(^/^|0f*('5i)l^/^) / dT,HmP-k?-<f)j^,\T,n{k)\\ (49) 

where T lh is the phase-space integral for the light hadrons and we only keep the leading 
dependence on momenta of the light hadrons. The formula is only for close to 4m^ and 
the predicted differential width is measurable in experiment. In the calculation we used 

(J/*|^Vx|0)(0|xV^|J/xl/) = {jmoi"'{'S,)\J/^>)e\e^)\ (50) 

In the equation e is the polarization vector of J/^P, the matrix element (J/\I'|Of^*('^S'i)| J/\E') 
is defined in and the average over the spin is implied in the matrix element. 

In this section we mainly have considered emission of gluons by a cc pair, and these 
gluons are transmitted into light hadrons. Because the kinematic region considered here, 
these gluons must be soft. That is why we have defined the nonperturbative object Tiuik) 
which includes gl- However, emission of hard gluons can happen, this can only happen as 
exchange of hard gluons between c-, c- quark and gluons. The effect of the exchange can 
be calculated with perturbative QCD, it will result in that the coefficient in Eq.(17) will be 
modified. For example, Z in Eq.(48) will be smaller than 1 by neglecting the effect at m~^. 
For exchange of soft gluons, the effect may be factorized into the NRQCD matrix element 
and Tlh- 

4. Summary and Discussion 
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In this work we have studied the exclusive J/\l/-decay into a lepton pair combined with 
hght hadrons, in which the total energy of the light hadrons is much smaller than Mj/$. 
The light hadrons are formed from soft gluons emitted by the c-flavored quarks. In the limit 
mc — >^ oo the nonperturbative effects related to J/^ and to the light hadrons can be sepa- 
rated and they are represented by a NRQCD matrix element and by a correlator of electric 
chromofields. Corrections to the hmit and those due to hard gluons may systematically 
added. We emphasize that our results are obtained without perturbative theory and are 
gauge invariant, although we firstly have derived them in the axial gauge with perturbative 
theory. But we have also derived them in an arbitrary gauge without using perturbative 
theory. 

In our approach an expansion in Vc, the velocity of the c-quark inside J/^ in its rest- 
frame, is used. We have only taken the leading order contributions at v^^. At this order, 
J/^ can be considered as a bound state of the cc pair in color singlet. The corrections to 
the leading-order results may also added. However, problems arise at order of f^. At this 
order J /'^ has a component in which the cc pair is in color octet and the component is a 
bound state of the cc pair with soft gluons. It is unclear how to add corrections from this 
component at order of v^. It deserves a further study of these problems to understand a 
bound state of many dynamical freedoms of QCD. 

Our approach applies for a class of decays. The results given here can be directly used 
for the T decay and can be generalized to B* ^ £u + light hadrons, where the light hadrons 
also have a total energy which is small. The same correlator appears in these decays to 
characterize the transition of soft gluons into light hadrons. In that sense the correlator is 
universal. Our results do not apply for 3^ because it is a spin-0 meson. However, with 
the approach presented here one may obtain corresponding predictions for Be- Our results 
do not apply for decays of excited quaukonium states like \1'' + £~ + tt + tt, in 

this decay after emission of soft gluons, which forms the light hadrons, the cc pair can be 
formed into J /"if, then the J/^' decays into a lepton pair. This type of contributions is not 
included in our results. But this type of contributions may be excluded by experimental 
cuts, through requiring the two leptons are not exactly back-to-back in the J/"^ rest-frame. 
With this requirement our results may apply. 

In this work we are unable to give any numerical result. In our prediction, the nonper- 
turbative effect related to J/^ is known, e.g., from the leptonic decay of J/^, while the 
nonperturbative effect related to the light hadrons is unknown. We currently develope a 
model for the correlator and will make an attempt to give numerical predictions. 

Acknowledgment: The author would like to thank Prof. Y.Q. Chen, Y.B. Dai, C. Liu, 
D. Soper and Z.X. Zhang for interesting discussions. This work is supported by National 
Science Foundation of P.R. China and by the Hundred Yonng Scientist Program of Sinica 
Academia of P.R.China. 



12 



REFERENCES 

[1] Z. G. Zhao, BES CoUabration, private communication 

[2] G.T Bodwin, E. Braaten, and, G.P Lepage, Phys. Rev. D51 (1995) 1125 

[3] N. Isgur and M.B. Wise, Phys. Lett. B232 (1989) 113, ibid B237 (1990) 527 
E. Eichten and B. Hill, Phys. Lett. B234 (1990) 511 
B. Grinstein, Nucl. Phys. B339 (1990) 253 
H. Georgi, Phys. Lett. B240 (1990) 447 

[4] J.C. Collins, L. Frankfurt and M. Strikman, Phys. Rev. D56 (1997) 2982 
J.C. Collins and A. Freund, Phys. Rev. D59 (1999) 074009 

[5] M.E. Peskin, Nucl. Phys. B156 (1979) 365, 

G. Bhanot and M.E. Peskin, Nucl. Phys. B156 (1979) 391 

[6] T.M. Yan, Phys. Rev. D22 (1980) 1652 

[7] M. Voloshin and V. Zakharov, Phys. Rev. Lett. Vol. 45 (1980) 688 
[8] O. Nachtmann, Ann. Phys. Vol. 209 (1991) 436 



13 



